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\  I.  Introduction. 

■V 

Friction  and/or  shear  have  been  suggested  by  many  people  as 
possible  sources  of  ignition  in  solid  explosives.  Although  these  two 
mechanisms  are  usually  considered  to  be  distinct,  on  a  microscopic 
scale  they  are  essentially  the  same.  On  a  macroscopic  scale,  friction 
involves  sliding  interfaces  and  is  described  in  terms  of  the  coeffi¬ 
cient  of  friction  and  the  stress  normal  to  the  interface.  The  rate  of 
heat  generation  at  the  surface  is  the  product  of  the  sliding  velocity 
times  the  coefficient  of  friction  times  the  normal  stress.  On  a  micro¬ 
scopic  level,  friction  is  caused  by  surface  irregularities  (asperities) 
which  deform  as  the  surfaces  slide.  The  frictional  heating  is  caused 
by  viscoplastic  work  on  the  asperities.  In  this  paper,  I  will  consider 
the  conditions  required  for  ignition  as  the  result  of  shear  defor¬ 
mation  and  viscoplastic  heating. 


$ 


A  cursory  examination  of  this  question  leads  one  to  the  con¬ 
clusion  that  the  deformation  must  be  localized  in  some  way  if  ignition 
is  to  occur.  The  plastic  work  per  unit  mass  resulting  from  the  defor¬ 
mation  of  a  material  is  given  by  the  following  relation  [1]: 
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deviatoric  strain,  e^  is  the  final  value  of  the  deviatoric  strain  and  p 

is  the  density.  Consider  the  homogeneous  deformation  of  a  cylinder  of 
explosive  so  that  the  length  decreases  and  the  radius  increases.  Assume 


Then  e^  =  eQ  = 


-1 


TV 


that  the  Tresca  yield  condition  holds. 

'r  -e  T  V  anG  °r  "  “0 
explosive,  and  the  usual  notation  of  cylindrical  coordinates  has  been 
used.  The  change  in  temperature  of  the  explosive  is  given  by 


S_  *  Sq  =  S^,  and  -  Sfl  =  Y,  where  Y  is  the  yield  strength  of  the 
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where  C  is  the  heat  capacity,  ZQ  is  the  initial  length,  and  z^  is  the 
final  length.  For  composition  B,  Y  is  about  0.07  GPa  and  C  is  about 

1.25  joules/°C.  To  obtain  a  100°C  change  in  temperature  requires 
zq/z  obtain  a  value  of  approximately  28.  Clearly  homogeneous  defor¬ 
mation  can  result  in  initiation  only  in  extreme  circumstances. 


There  are  several  circumstances  which  may  lead  to  very  high 
shears  in  localized  regions.  Some  of  them  are  the  following: 

1.  Sliding  friction.  As  explained  above,  on  a  microscopic 
scale  sliding  friction  involves  shear  deformation,  and  very  high 
strains  may  occur  in  narrow  layers.  In  particular,  friction  between 
a  rotating  shell  and  the  explosive  fill  has  been  suggested  as  a  cause 
of  premature  ignition  in  artillery  shells  [2],  It  has  also  been  sug¬ 
gested  as  a  mechanism  for  the  propagation  of  reaction  during  nondeto- 
native  explosions  in  confined  explosive  charges  [3] . 


2.  The  collapse  of  spherical  cavities.  During  ductile 
cavity  collapse,  a  thin  layer  of  material  at  the  boundary  of  the 
cavity  sees  very  high  strains.  This  situation  has  been  considered 

by  Carroll  and  Holt  [4].  Their  analysis  indicates  that,  in  the  absence 
of  melting,  extremely  high  temperatures  car  be  obtained  in  small  vol¬ 
umes  around  a  collapsed  cavity.  This  has  been  suggested  as  an  igni¬ 
tion  mechanism  during  shock  wave  initiation  of  explosives. 

3.  Adiabatic  shear  bands.  Adiabatic  shear  bands  in  metals 
have  been  discussed  by  Recht  [5].  A  shear  band  is  likely  to  form  in 
a  deforming  material  when  thermal  softening  exceeds  work  hardening. 

In  this  case,  deformation  in  any  plane  causes  the  material  in  that 
plane  to  weaken  and  tends  to  concentrate  further  deformation  in  the 
same  plane.  In  metals,  thin  bands  of  material  can  be  heated  to  the 
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melting  point.  Winter  and  Field  [6]  have  applied  Recht’s  analysis  to 
explosives  and  concluded  that  explosives  should  be  highly  susceptible 
to  shear  banding.  Afanasev  [7]  and  co-workers  have  observed  them  in 
nonexplosive  organics. 

4.  Extrusion  of  explosive^ into  cracks.  If  a  crack  develops 
in  the  metal  confinement  around  an  explosive  while  the  explosive  is 
under  pressure,  ignition  could  result  when  the  explosive  extrudes  into 
the  crack.  This  has  been  suggested  as  a  possible  ignition  mechanism 
when  confined  explosives  are  impacted  by  projectiles  at  velocities 
too  low  to  cause  shock  wave  initiation  [3] . 

In  all  of  these  circumstances,  melting  tends  to  limit  the 
temperature  which  can  be  achieved.  Once  the  melting  point  has  been 
reached,  any  further  heating  must  be  strictly  by  viscous  processes. 

The  melting  point  increases  with  pressure,  and  melting  is  a  more  impor¬ 
tant  consideration  at  low  pressure  than  at  high  pressure.  From  data 
in  the  International  Critical  Tables,  the  rate  of  increassoof  melting 
point  with  pressure  for  many  organic  materials  is  20  to  30°  per  ^i lo¬ 
bar.  For  instance,  for  TNT  the  melting  point  increases  about  25°/kilo- 
bar  at  atmospheric  pressure.  The  rate  of  increase  decreases  as  pres¬ 
sure  increases.  To  obtain  a  crude  estimate  of  the  effect  of  melting 
on  hot  spot  ignition,  I  assumed  £hat  the  melting  point  of  several 
explosives  increased  at  about  20°/kilobar  and  computed  the  time  to 
explosion  at  the  melting  point  as  a  function  of  pressure.  Kinetic 
constants  were  taken  from  Rogers  [8] .  The  results  are  shown  in  Figure 
1.  For  reactions  which  occur  in  1  microsecond  or  less,  which  is 
typical  of  the  shock  initiation  regime,  melting  is  important  in  pre¬ 
venting  initiation  at  pressures  below  1.0  GPa  in  PETN,  RDX,  or  HMX  and 
for  pressures  below  3.0  GPa  in  TNT.  For  reactions  which  occur  in  a 
millisecond  time  frame,  which  might  be  typical  of  artillery  setback 
conditions,  melting  is  significant  for  pressures  below  0.4  GPa  in  RDX 
and  1.5  GPa  in  TNT.  In  this  paper  I  will  investigate  the  conditions 
under  which  viscoplastic  heating  can  lead  to  ignition  when  melting 
is  considered. 

II.  Ductile  Cavity  Collapse. 

Carroll  and  Holt  treated  spherical  cavity  collapse  in  an 
elastic-plastic  fashion  [4]  and  in  an  elastic-viscoplastic  fashion  [9]. 
Application  of  their  equations  shows  that,  in  the  absence  of  melting, 
extremely  high  temperatures  are  produced  in  a  thin  layer  around  the 
collapsing  cavity.  Figure  2  shows  the  temperature  as  a  function  of 
applied  pressure  for  a  Comp-B  charge  where  the  initial  porosity  (actual 
volume/fully  compacted  volume)  is  1.03.  Extremely  high  temperatures 
are  predicted  in  a  very  thin  layer  around  the  collapsing  cavity.  At  a 
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pressure  of  g.O  GPa,  the  surface  of  the  cavity  reaches  a  temperature 
of  about  600°C.  However,  only  a  very  thin  layer  is  significantly 
heated.  Figure  2  also  shows  the  temperature  of  the  isotherms  enclosing 
the  hottest  0,01%  and  the  hottest  1%  of  the  material.  These  isotherms 
are  much  cooler. 


The  analysis  of  Carroll  and  Holt  does  not  include  thermal 
softening,  melting,  or  heat  conduction.  Since  the  heated  region  is 
close  to  the  surface  of  the  cavity,  where  the  pressure  is  low,  one 
can  suspect  that  the  inclusion  of  these  factors  might  alter  the 
results  significantly.  We  have  modified  Carroll  and  Holt's  formu¬ 
lation  slightly  to  include  these  effects  and  also  added  a  viscosity 
which  varies  with  temperature  and  pressure.  For  simplicity,  we  use  a 
rigid-plastic  model  (i.e.,  we  assume  that  the  elastic  strains  are  negli¬ 
gible).  The  principal  assumptions  of  the  Carroll  and  Holt  model  are 
the  following: 


symmetric. 


1.  The  flow  in  the  vicinity  of  the  cavity  is  spherically 


2.  The  material  is  incompressible. 

3.  The  following  constitutive  relation  applies: 


_L  s 

2G  bij 


(1  -  Y/vST)  S. 


where  is  a  component  of  the  deviator  strain  rate  tensor,  S... 

is  a  component  of  the  deviator  stress  tensor,  G  is  the  shear  modulus, 
y  is  the  viscosity,  Y  is  the  yield  strength,  and  I  is  the  second 
invariant  of  the  deviator  stress.  For  the  spherically  symmetric, 
rigid-plastic  case  considered  here,  this  reduces  to 


Sr  =  2yer  +  3  Y 


se  s  2m60  3  Y  ’ 

where  Sr,  Sg,  er>  eg  refer  to  the  principle  deviator  stresses  and 

principle  deviator  strain  rates,  and  the  usual  notacior  for  spherical 
coordinates  has  been  used. 

Carroll  and  Holt  considered  a  hollow  sphere  with  outer 
radius  b,  inner  radius  a,  and  porosity  a,  defined  as  the  ratio  of  the 


u.*  v 


46 


k  i 


FREY 


3  3  3 

actual  volume  to  the  fully  compacted  volume,  where  a  =  b  /(b  -a  ), 
The  equation  of  motion  for  this  system  is 

3a 


‘r  2 

T"  +  r 
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°e> 


pr  , 


where  ar  and  cQ  are  the  principle  stresses,  r  is  the  radius,  and  dots 

imply  differentiation  with  respect  to  time.  Following  Carroll  and 
Holt,  we  note  that 

ar  -  a9  =  Sr  ■  Se  =  u  er  +  Y  * 
and  apply  the  boundary  conditions, 

a  =  0  at  r  =  a 
r 

a  =  -P  at  r  =  b, 
r  * 

where  P  is  the  applied  pressure. 


We  also  note  that 
2  B 


e  =  +  3 

r  3  r 

where  B(t)  =  a3  -  a3  = 


b3  -  b3  =  a  3  (a  ~a)/(a  -1) 
O  O  v  0  '  v  o 


and  aQ,  bQ,  and  aQ  are  the  initial  values  of  a,  b,  and  a.  The  equation 

of  motion  can  then  be  integrated  with  respect  to  the  radius  and  the 
result  expressed  as  the  following  differential  equation  for  a: 
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The  rate  of  plastic  work  at  any  point  in  the  flow  may  be  computed  from 
equation  (1),  giving  the  following  result: 


Wp  CD  -  I  V 


(3) 
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Where  is  the  rate  of  plastic  work  per  unit  volume.  We  have  :'.nte- 

grated  equation  (2)  numerically.  Simultaneously,  we  integrate  the 
heat  equation  for  the  material  around  the  cavity. 


3T 

at 


_K  3^T  2  3T  *£ 

pC  3r2  +  r  8r  pC 


(4) 


where  T  is  temperature,  t  is  time,  K  is  heat  conductivity,  p  is  density, 
and  C  is  heat  capacity.  This  equation  was  also  integrated  using  finite 
differences  and  a  computational  mesh  which  concentrated  most  of  the 
points  in  the  region  close *to  the  surface  of  the  cavity.  At  each  time 


step, equation  (2)  was  integrated  to  obtain  a,  a,  a,  a,  and  b.  Then  W 


could  be  determined  and  equation  (4)  integrated  to  determine  T(r)  at 
the  next  time  step. 


Good  data  on  the  temperature  dependence  of  the  yield  strength 
and  on  the  temperature  and  pressure  dependence  of  the  viscosity  of 
common  explosives  are  not  available.  We  assumed  that  the  melting  point 
increased  with  pressure  at  the  rate  of  200/GPa  and  that  the  yield 
strength  was  zero  above  the  meltigg  point  and  decreased  to  zero  in  a 
linear  fashion  over  a  range  of  30°C  below  the  melting  point.  We 
assumed  that  the  viscosity  varied  in  the  following  way  with  temperature 
and  pressure: 


where  yQ  is  the  viscosity  at  temperature  T  ,  E  is  an  experimentally 
determined  activation  energy  for  viscous  flow,  and  Pq  is  an  experi¬ 
mentally  determined  constant.  Pressure  and  temperature  relations  of 
this  type  are  discussed  by  Frenkel  [10]  and  have  some  theoretical 
support.  The  pressure  dependence  shown  in  equation  (5)  agrees  well 
with  data  from  the  International  Critical  Tables.  Some  data  is 
available  for  TNT,  and  we  will  use  this  material  in  all  of  the  calcu¬ 
lations  reported  here.  We  determined  an  approximate  value  for  E, 
3,880u  K,  on  the  basis  of  two  data  points  [11].  We  estimated  PQ  to 

be  0.165  GPa,  a  value  which  is  typical  of  organic  liquids  and  fits 
the  experimental  data  reported  in  the  International  Critical  Tables 
for  nitrobenzene.  Other  required  parameters  are  shown  in  Table  I, 
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TABLE  I  -  MATERIAL  PROPERTIES  USED  IN  CALCULATIONS 


Parameter 

Y 

v 

o 

C 

K 

P 

E 


_ Value _ 

0.07  GPa 

1.39  x  10_2Kg/ (msec)  at  85°C35 
1.29  x  105  joule/(Kg°C)35 
0.262  joule/ (m°C  sec)33 
1.64  x  103  Kg/m3  35 
3, 880°K 


P  0.165  GPa 

o 

Figure  3  shows  the  results  of  the  calculation  for  a  case 
where  the  pressure  increased  from  zero  to  2.0  GPa  over  a  period  of 
20  nanoseconds  (this  simulates  a  shock  wave  in  a  porous  medium) .  The 
peak  temperature  never  exceeds  the  normal  (one  atmosphere)  melting 
point.  In  this  situation,  the  deformation  is  concentrated  in  the 
region  close  to  cavity  where  the  pressure  and  melting  point  are  low. 
Viscosity  is  not  sufficient  to  elevate  the  temperature  significantly 
above  the  melting  point.  Therefore,  ignition  by  this  mechanism  is 
unlikely  in  TNT,  Although  we  have  not  performed  calculations  on 
other  explosives,  we  can  infer  that  they  would  not  be  heated  much 
above  the  normal  melting  point  either.  Since  RDX  has  a  rather  long 
thermal  explosion  time  at  the  melting  point,  it  is  probably  not 
susceptible  to  ignition  by  this  mechanism  in  a  millisecond  time  scale. 


Ill,  Friction  and/or  adiabatic  shear  bands. 


We  have  already  commented  on  the  similarity  of  friction  and 
adiabatic  shear  bands  at  the  microscopic  level.  In  each  case,  large 
deformations  are  concentrated  in  a  narrow  region.  In  a  shear  band,  the 
concentration  occurs  because  thermal  softening  exceeds  work  hardening, 
as  explained  above.  At  a  sliding  interace,  the  concentration  occurs 
because  the  interface  region  is  weaker  than  the  rest  of  the  material. 

We  have  calculated  the  time  evolution  of  temperature  in  a  shear  band 
or  at  a  sliding  interface  when  melting  occurs.  Our  intention  was  to 
perform  the  simplest  possible  calculation  while  retaining  the  essential 
features  of  the  real  world. 
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In  the  model,  we  impose  a  velocity  gradient  across  an  arbi¬ 
trary  thickness  of  explosive.  The  magnitude  of  this  thickness  does 
hot  influence  the  results  of  the  calculation.  We  arbitrarily  induce  a 
shear  band  by  specifying  the  initial  temperature  in  a  thin  layer  to 
be  high  enough  so  that  thermal  softening  occurs.  The  time  evolution 
of  the  system  is  followed  by  making  use  of  the  following  assumptions: 


1.  The  shear  stress  is  constant  in  all  planes  parallel  to 
the  shear  band  at  any  instant;  i.e.. 


u  (x) 


dv  +  cf(x)  =  f(t) 

17 


(6) 


where  y  is  viscosity,  x  is  distance  perpendicular  to  the  shear  band 
(perpendicular  to  the  direction  of  flow),"v  is  velocity,  a  is  shear 
strength,  and  f  is  the  shear  stress  applied  at  the  shear  band,  f  is 
independent  of  x  but  may  vary  with  time.  This  assumption  is  a  quasi¬ 
steady  state  assumption.  It  requires  that  the  flow  at  any  instant 
must  have  the  form  it  would  have  at  steady  state  given  the  temperature 
distribution  which  exists  at  that  time,  but  it  permits  the  flow  to 
change  as  the  temperature  distribution  changes. 


2.  The  total  shear  velocity  (the  integral  of  dv/dx  with 
respect  to  x)  is  independent  of  time.  This  condition  and  the  previous 
one  determine  f  via  the  following  equation: 

/k  dx  ■  /lif- dx  ■  vt  •  (,) 

where  vt  is  the  total  shear  velocity.  At  any  time  if  a(x)  and  y(x)  are 

known,  this  equation  can  be  used  to  determine  f,  and  equation  (6)  can 
be  used  to  compute  dv/dx  as  a  function  of  x.  As  the  shear  band 
develops,  the  velocity  gradient  will  become  concentrated  in  a  narrow 
region  and  dv/dx  will  be  ■'.ero  over  most  of  the  space.  In  this  situa¬ 
tion,  the  integrals  in  equation  (7)  must  be  evaluated  only  over  the 
region  where  dv/dx  is  positive.  This  is  accomplished  by  choosing  the 
limits  on  the  second  integral,  which  is  evaluated  numerically,  so  that 
the  shear  stress  f  is  minimized. 


3.  The  solid  is  rigid-perfectly  plastic  with  a  shear 
strength  which  decreases  to  zero  at  the  melting  point.  As  in  the 
previous  discussion,  the  shear  strength  is  assumed  to  decrease  linearly 
to  zero  in  a  30°C  range  of  temperature  below  the  melting  point. 

4.  The  viscosity  varies  with  temperature  and  pressure  in  the 
manner  described  by  equation  (5) . 
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5.  The  rate  of  viscoplastic  work  per  unit  volume  is  given 
by  the  fol lowing  expression: 


6.  The  heat  generated  per  unit  mass  by  chemical  decom¬ 
position  can  be  computed  on  the  basis  of  simple  Arrhenius  kinetics, 
as  follows: 


Q  =  QA  exp  (-Ea/T) 

where  Q  is  the  heat  of  reaction  per  unit  mass,  A  is  the  frequency 
•  factor,  and  E  is  the  activation  energy. 

a 

7.  One  dimensional  heat  equation  applies;  i.e., 


?  *  w 

3T  JC  jTT  Q  JL 

*t  "  pC  3x2  +  C  +  pC  * 


(8) 


Equation  (7)  was  integrated  numerically  using  finite 
differences.  At  each  time-step,  a(x)  and  p(x)  were  determined  from 
the  known  temperature  distribution  and  dv/dx  was  determined  from 
equations  (6)  and  (7).  Melting  is  handled  in  one  of  two  ways.  In 
one  procedure,  all  of  the  material  is  considered  to  be  a  single  phase, 
but  the  shear  strength  goes  to  zero  when  the  temperature  exceeds  the 
melting  point.  This  procedure  neglects  the  heat  of  fusion.  In  the 
other  procedure,  two  phases  and  a  phase  boundary  are  included  in  the 
calculation.  The  velocity  of  the  phase  boundary  (the  rate  of  growth 
of  the  melt  layer)  is  detennined  from  an  energy  balance  at  the 
boundary: 

3T  9T 

VQf =  Ke  93T  *  Ks  inr  * 


where  is  the  velocity  of  the  boundary,  is  the  heat  of  fusion,  and 

subscripts  e  and  s  refer  to  the  liquid  and  solid  phases.  When  this 
procedure  is  used,  a  dynamic  mesh  is  used  in  both  phases  to  avoid  the 
necessity  of  rezoning  as  the  boundary  moves.  This  is  accomplished  by 
introducing  a  variable,  y,  defined  as  follows: 


where  xbQ  is  the  initial  position  of  the  boundary  and  xb  is  the  position 

of  the  boundary  at  any  time,  x  and  its  derivatives  ai‘e  expressed  in 
terms  of  Y  and  substituted  into  equation  (8).  The  resulting  equation 


is  solved  by  finite  differences.  When  expressed  in  terms  of  y,  the 
computational  grid  automatically  expands  or  contracts  to  adjust  to 
the  moving  boundary.  No  significant  differences  were  noted  in  the 
results  of  the  two  types  of  calculation.  Surprisingly,  the  second  pro¬ 
cedure  gave  slightly  higher  temperatures,  but  this  may  be  due  to 
numerical  error.  All  of  the  results  reported  here  were  accomplished 
with  the  second  procedure. 

Calculations  were  performed  using  material  properties  for 
TNT  as  listed  in  Table  I.  In  addition,  the  following  values  for 
Q,  A,  and  E  were  taken  from  Rogers  [8] : 

Q  =  300  Cal/ (°C  g) 

A  --  2.51X1011  sec"1 

Ea=  34.4  KCal/mole 

Since  the  model  is  only  an  approximation  of  reality  and  since  many 
parameters  are  not  known  precisely,  the  results  should  be  considered 
in  a  qualitative  rather  than  a  quantitative  sense.  Figure  4  shows 
how  the  computed  temperature  at  the  center  of  the  shear  band  varies 
with  time  for  a  case  where  the  shear  velocity  was  0.2  km/sec  and  the 
pressure  was  1.03  GPa.  The  temperature  asymptotically  approaches  a 
maximum  value.  The  time  required  for  significant  heating  to  occur 
is  very  short.  The  initial  width  of  the  heated  region,  which  triggers 
shear  band  formation,  affects  the  time  history  at  the  shear  band  but 
does  not  affect  the  maximum  temperature.  Figure  5  shows  the  time 
history  for  several  calculations  with  the  same  pressure  and  shear 
velocity  but  different  initial  widths.  There  is  an  optimum  initial 
width  which  gives  the  most  rapid  heating  rate.  With  smaller  initial 
widths,  the  temperature  in  the  shear  band  may  decrease  for  a  while 
before  it  begins  to  increase.  For  larger  initial  widths,  the  tem¬ 
perature  increases  more  slowly.  In  the  absence  of  reaction,  the 
maximum  attainable  temperature  is  a  strong  function  of  pressure. 

Figure  6  shows  how  the  maximum  temperature  varies  with  shear  velocity 
and  pressure.  At  low  pressures,  melting  limits  the  attainable  tem¬ 
perature  even  for  very  high  shear  velocities.  High  pressure  increases 
both  the  melting  point  and  the  viscosity  and  permits  much  higher  tem¬ 
peratures  to  be  attained  at  lower  shear  velocities. 

On  the  basis  of  these  calculations,  it  appears  very  unlikely 
that  friction  couxd  be  responsible  for  a  premature  ignition  in  a  TNT 
shell  during  gun  launch.  The  maximum  pressure  during  gun  launch  is 
less  that  0.2  GPa,  and  the  sliding  velocity  between  the  rotating  shell 
and  the  explosive  is  less  than  30  m/sec.  Under  these  conditions,  the 
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temperature  would  not  rise  significantly  above  the  melting  point  (of 
course,  the  presence  of  high  melting  point  grit  would  change  the  pic¬ 
ture  drastically),  hfe  have  not  done  any  calculations  on  RDX  because 
we  do  not  have  data  for  its  viscosity.  However,  we  can  make  a  rough 
estimate  of  what  would  happen  in  RDX  by  using  the  TNT  parameters  but 
adjusting  the  melting  point  and  T  (from  equation  (5)  for  the  vis¬ 
cosity)  upwards  to  account  for  the  difference  in  melting  point  between 
TNT  and  RDX.  If  we  do  this,  we  obtain  the  same  curves  which  we  ob¬ 
tained  for  TNT  but  with  all  of  the  temperatures  higher  by  124  C  (the 
difference  in  melting  point  between  RDX  and  TNT).  On  this  basis,  it 
appears  that  RDX  should  also  be  insensitive  to  friction  under  setback 
conditions.  These  conclusions  should  be  checked  by  experiment,  and  a 
program  to  do  so  is  currently  in  progress  at  BRL.  At  present,  the 
experiments  support  the  theoretical  conclusi  u. 

Under  shock  initiation  conditions ,  shear  bands  would  be 
expected  to  form  in  the  vicinity  of  collapsing  cavities.  Figure  7 
shows  in  a  schematic  fashion  how  this  might  happen.  In  this  situation, 
both  the  pressure  and  the  shear  velocity  could  be  much  higher  than  in 
the  setback  situation.  An  estimate  of  the  maximum  shear  velocity  can 
be  obtained  by  taking  the  difference  between  the  free  surface  velocity 
and  the  particle  velocity  behind  the  shock.  This  difference  is  approx¬ 
imately  equal  to  the  particle  velocity.  A  1.03  GPa  shock  in  TNT  has  a 
particle  velocity  of  about  0.2  km/sec.  Figure  8  shows  how  the  computed 
shear  band  temperature  varies  with  time  in  this  case.  Thermal  explosion 
occurs  after  about  SO  nanoseconds.  However,  if  the  shear  stops  before 
thermal  explosion  occurs,  the  temperature  in  the  shear  band  decreases 
rapidly  due  to  the  steep  thermal  gradients,  and  reaction  does  not  occur. 
This  is  also  shown  in  Figure  8.  An  estimate  of  the  maximum  slip  dis¬ 
tance  at  a  shear  band  of  this  type  is  the  dimension  of  the  cavity  in 
the  direction  of  shock  motion.  The  duration  of  the  shear  is  this  dis¬ 
tance  divided  by  the  shear  velocity.  Consequently,  we  can  estimate, 
ver>  approximately,  the  critical  cavity  size  required  for  initiation  by 
this  mechanism  at  any  shock  pressure.  Table  II  shows  these  estimates 
for  cavities  in  TNT. 

TABLE  II  -  CRITICAL  CAVITY  SIZES  FOR  SHOCK  IGNITION  IN  TNT, 

COMPUTED  ON  THE  BASIS  OF  SHEAR  BAND  IGNITION 


Shock 
Pressure 
(K  bar) 

Particle 

Velocity 

(Km/sec) 

Time  to 

Explosion 

(Microsecond) 

Critical 

Cavity 

Size 

4,9 

0.1 

4700 

4.7cm 

7,5 

0.15 

1.75 

0. 26mm 

10.3 

0.20 

0.05 

0. 1mm 
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Liddiard’s  data  [12]  showed  a  critical  pressure  of  about  0.7  GPa  for 
the  ignition  of  TNT  by  long  duration  shocks.  Since  the  maximum  cavity 
size  in  a  TNT  charge  of  reasonable  quality  would  be  about  0.1  to  0.3 
tan,  the  computed  numbers  are  in  reasonable  agreement  with  experiment. 

IV.  Summary  and  Conclusions. 

The  role  of  shear  in  the  initiation  of  reaction  in  high  explo¬ 
sive  has  been  examined  theoretically.  The  calculations  show  that,  for 
initiation  by  shear  to  occur  under  realistic  conditions,  some  local¬ 
ization  of  the  deformation  must  occur.  When  localized  shear  occurs, 
pressure  is  a  critical  parameter  which  controls  the  occurrence  of 
ignition.  High  pressure  during  shear  deformation  enhances  the  chances 
for  ignition  by  raising  both  the  melting  point  (thereby  increasing  the 
heating  due  to  viscoplastic  work)  and  the  viscosity  of  Che  melt  (hence, 
increasing  viscous  heating) .  This  pressure  effect  is  extremely  impor¬ 
tant  in  causing  an  explosive  to  be  sensitized  to  shear  deformation, 
and  should  be  kept  in  mind  when  one  assesses  explosive  hazards.  For  a 
low  melting  explosive,  such  as  TNT,  pressures  in  excess  of  0.5  GPa  are 
required  for  ignition  at  any  shear  velocity  likely  to  occur.  Fox 
higher  melting  explosives,  lower  pressures  can  lead  to  initiation. 

The  calculations  show  that  typical  pressures  encountered 
during  gun  launch  are  unlikely  to  cause  ignition  of  either  TNT  or  RDX 
by  this  mechanism.  However,  when  ignitions  occur  as  a  result  of  other 
causes,  the  resulting  higher  pressures  make  this  mechanism  a  serious 
potential  source  of  secondary  ignitions  which  lead  to  violent  explo¬ 
sions.  This  mechanism  is  also  considered  to  be  an  important  pathway  by 
which  interround  propagation  occurs,  when  munitions  within  a  storage 
array  are  subjected  to  multiple  fragment  impacts  and  severe  crushing. 
(See  reference  3.) 

Ductile  cavity  collapse  does  not  lead  to  conditions  which 
cause  ignition.  During  ductile  cavity  collapse,  the  strain  is  con¬ 
centrated  in  regions  where  the  pressure  and  viscosity  are  low,  and  the 
temperature  which  may  be  obtained  is  limited  by  the  melting  point. 
Brittle  cavity  collapse,  such  a?  might  be  expected  under  shock  loading 
conditions,  can  lead  to  the  formation  of  shear  bands  and  high  tem¬ 
peratures  can  be  attained.  This  is  also  expected  to  occur  when  the 
explosive  is  extruded  into  cracks  in  the  confining  steel  case  while 
the  munition  is  under  pressure.  This  latter  sequence  of  events  is 
viewed  as  a  very  probable  source  of  ignition  when  the  casing  fails 
under  launch  or  when  subjected  to  crushing. 


54 


VQ/* 


i  X 


Figure  1  -  Thermal  Explosion  Time  at 
the  Melting  Point  vs  Pressure. 


Figure  3  -  Temperature  vs  Porosity 
Durin;  Cavity  Collapse  When  Melting 
and  Heat  Conduction  Occur. 
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Figure  2  -  Temperature  vs  Pressure 
After  Cavity  Collapse  (Carroll  5 
Holt's  Model). 
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Figure  4  -  Temperature  vs  Time  in  a 
Developing  Shear  Band. 
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Figure  S  -  The  Time  History  of  a  Shear 
Sand  as  a  Function  of  the  Width  of 
the  Initial  Heated  Region. 


Figure  6  -  The  Effect  of  Shear  Velocity 
and  Pressure  on  Peak  Temperatures  in  a 
Shear  Band. 
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Figure  7  -  Schematic  Drawing 
Shear  Bands  Could  Be  Formed 
Explosive. 
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Figure  3  -  Temperature  vs  Time  in  a  Shear 
Band. 
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